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We report the first measurement of tlie excitation spectrum and the static structure factor of a Bose- 
Einstein condensate. The excitation spectrum displays a linear phonon regime, as well as a parabolic 
single-particle regime. The linear regime provides an upper limit for the superfluid critical velocity, 
by the Landau criterion. The excitation spectrum agrees well with the Bogoliubov spectrum, in 
the local density approximation. This agreement continues even for excitations close to the long- 
wavelength limit of the region of applicability of the approximation. Feynman's relation between 
the excitation spectrum and the static structure factor is verified, within an overall constant. 



The excitation spectrum of superfluid ^He gives im- 
portant insights into this superfluid The excitation 
spectrum places an upper bound on the superfluid criti- 
cal velocity. Furthermore, it indicates the types of exci- 
tations which occur in the superfluid, and reflects the su- 
perfluid's density correlations. The excitation spectrum 
of a Bose-Einstein condensate should give similar insight 
into this system. 

The excitation spectrum is the relation uj{k), giving 
the energy Tiuj(k) of each excitation, as a function of its 
wave vector k. For excitations with wavelengths 27r/fc 
which are comparable to the radius of the condensate, the 
spectrum is characterized by discrete oscillatory modes, 
which are strongly dependent on the shape of the conden- 
sate. These modes have been measured previously 

For wavelengths much shorter than the radius of the 
condensate in the direction of fc, the excitation spectrum 
becomes an essentially continuous function of fc, which 
characterizes the intrinsic bulk properties of the conden- 
sate In this work, "excitation spectrum" refers to 
this bulk regime. We report the flrst measurement of 
the excitation spectrum. 

Previously, for condensates with various values of the 
chemical potential /i, isolated points on the excitation 
spectrum were measured J^, f^. These points, 

and their dependence on /i, were found to be consistent 
with Bogoliubov theory. 

For a homogeneous condensate (constant density n), 
the zero-temperature excitation spectrum is expected to 
be of the Bogoliubov form which is valid as long as 
the wavelength of the excitation is much greater than the 
s-wave scattering length a. The Bogoliubov excitation 
spectrum is given by 
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The speed of sound c is given by gn/m, where the 
constant g is given by A-Kfi^a/m and m is the atomic 
mass. The chemical potential ^ = gn characterizes the 
strength of the atomic interactions. In the limit of small 
k, u){k) w ck. This linear part of the excitation spectrum 
corresponds to phonons. 



By definition, uj{k) is the average frequency of the dy- 
namic structure factor ^ S{k,Lu), which gives the re- 
sponse of the excitation process. Integrating S^kjUj) 
over uj gives the static structure factor S{k), which is 
the Fourier transform of the density correlation function, 
giving the magnitude of the density fluctuations / [J^/ in 
the fluid, at wavelength 2'K/k. In general, the /-sum 
rule ]^ yields Feynman's relation J7^/ 

S{k)^{hk^/2m)/uj{k) (2) 

Eq. ^ relates the strength of the resonance to its fre- 
quency. 

For a BEC in a parabolic trap, the density is inhomo- 
geneous, so does not precisely apply. However, we 
will see that the form of ijj{k) is very close to that of (|^). 
In the local density approximation ^ (LDA), the speed 
of sound, and therefore the excitation spectrum (^, can 
be considered to be deflned at each point r in the BEC, 
where c = gn{f)/m. This approximation is valid as 
long as the Thomas-Fermi radius of the condensate in 
the fc-direction is much larger than the wavelength of the 
excitation J^, ^l^j, JT^. For the parabolic density 
proflle of the condensate in the Thomas-Fermi regime, 
the excitation spectrum for the entire condensate in the 
LDA is 
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where ciaik) is given by ^^S'(fc)-2 - 1, and ^, ^ 
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where a = 2fi/[h'^k'^ /{2m)] and ^ = gn{0), where 
n(0) is the maximum density in the condensate. The 
value cid {k) is a weak, monotonically increasing function, 
which varies from Ceff = 32/(157r)-y/ fj,/m{= 0.68 -^//i/m) 
for small fc, to ciarge — •\/4/7\/ /i/m(= 0.76^/ /i/m) for 
large k. The value ciarge is just the speed of sound for a 
homogeneous condensate whose density is the density of 
the inhomogeneous condensate, averaged over each atom. 
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Since cid{k) is nearly constant, the excitation spectrum in 
the LDA (HI) is very similar in form to the homogeneous 
excitation spectrum (|lj), illustrating that the excitation 
spectrum of the trapped condensate should largely reflect 
the intrinsic properties of the homogeneous condensate. 
Thus, the functional form of aj(fc) reported here provides 
a test for the Bogoliubov theory for a homogeneous con- 
densate. 

Since small k corresponds to phonons, Ceff is the ef- 
fective speed of sound for an inhomogeneous condensate. 
The Bogoliubov spectrum in the LDA fulfills the require- 
ment of Bose statistics that the lowest lying excitations 
be phonons JT^/ . Phonons are collective excitations, each 
of which consists of a large number fT^I Nk of atoms 
moving with momentum hk, and Nk — 1 atoms moving 
with momentum -hk. This momentum distribution of a 
phonon was measured in Previously, a sound pulse, 
composed of a combination of phonons, was imaged in a 
time sequence of absorption images J7^. The observed 
sound velocity was in rough agreement with Ce/ / . 

For large fc, (||) is given by 
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where the first term is much larger than the second. This 
parabolic part of the spectrum corresponds to single- 
particle excitations, in which a single atom has a veloc- 
ity hk/m, which is much larger than the speed of sound. 
The second term in (||) is equal to {4:/7)fi/h and is inde- 
pendent of fc, reflecting the extra interaction energy ^ 
experienced by a moving atom JT^. 

The transition from collective excitations to single- 
particle excitations occurs for k on the order of 
where ^ is the healing length /p^. We take to 
be the solution of fc = ^/2mcid{k) /Ti, which is given by 

= \/2mc/Tl, and C « (Ce// -|- Clarge)/'^- 

An upper limit on the superfluid critical velocity can 
be extracted from the excitation spectrum, via the Lan- 
dau criterion, which states that the superfluid cannot 
flow with a speed greater than w/fc, for any excitation 
Lo{k) in the spectrum Therefore, the critical veloc- 
ity Vc cannot be more than the smallest value oi oj /k in 
the entire spectrum. According to Bogoliubov theory, 
the phonons have the smallest value oi uj/k, as seen in 
Fig. 2a, so Vc cannot be greater than the speed of sound 
Cef / . We emphasize that the critical velocity for vortex 
production is usually much lower than the critical veloc- 
ity given by the Landau criterion, so vortex production 
usually limits the speed of superfluid flow J7^/ , jl^. 

The condensate of *^Rb atoms in the 5si/2, F = 2, 
niF = 2 ground state is produced in a QUIC magnetic 
trap loaded by a double MOT system. The mag- 
netic trap contains 6x10^ atoms. After 22 seconds 
of evaporation, 1x10^ atoms remain, forming a nearly 
pure condensate, with a thermal fraction of 5% or less. 
The bias magnetic field is 2 G. The radial and axial 



trapping frequencies are 25 Hz and 220 Hz respectively, 
yielding radial and axial Thomas-Fermi radii of 3 
and i? = 28 fj,m, respectively. The chemical potential 
fi/h = 2.02 ± 0.09 kHz is determined ^ by the radial 
size of the condensate in a time-of-flight image. 

The excitation spectrum is measured by Bragg spec- 
troscopy J2^. Two Bragg beams A and B with 
approximately parallel polarization Jl^, separated by an 
angle 3° < 6* < 130°, illuminate the condensate for a 
time ts- The frequency of beam A is greater than the 
frequency of beam B by an amount lu determined by two 
acousto-optic modulators. If a photon is absorbed from 
A and emitted into i?, an excitation is produced with en- 
ergy LU and momentum k, where k = 2kp sin {0/2), and kp 
is the photon wave number. Here, we neglect the possi- 
bility that a single photon will excite multiple excitations, 
in contrast to the case of superfluid *He, in which mul- 
tiparticle excitations are an important contribution j^/ , 
^to S{k,u}). 

The wave vector k is precisely adjusted to be along the 
axis of the cigar-shaped condensate Jj^/ . To insure that 
the entire condensate is stimulated by the Bragg pulse, 
the length of the pulse is is chosen such that the spectral 
width of the pulse is roughly equal to the intrinsic width 
of the resonance. For this experiment, the broadening 
due to inhomogeneous density l^vid always dominates the 
Doppler broadening and is given by ;§/ 0.45 kHz for 
large k, and 0.3 ti;(fc)/(27r) in the phonon regime. Thus, 
ts is chosen to be roughly {2I^vid)~^ ■ For large fc, the 
resonance may be further broadened by s-wave scattering 

m 

The beams are detuned A = 6.5 GHz below the 5Si/2, 
5P3/21 F — i transition. The intensities I a 
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and Ib of each beam are adjusted to values between 0.1 
mW cm~^ and 1.1 mW cm~^ so that the number 
of excitations is 10% to 20% of the number of atoms in 
the condensate. For pulses of this strength, the chemical 
potential decreases by an average of only 12% during the 
pulse, so the Bragg scattering can be considered to be a 
small perturbation to the condensate. The time average 
of ji/h during the pulse is 1.91 ±0.09 kHz, which is taken 
to be the relevant value for the excitation spectrum. 

To measure a single point on the excitation spectrum 
Lu{k), k is fixed by 0, and to is varied. The resonant value 
of OJ is taken as uj{k). 

After the Bragg pulse, the atoms are allowed to expand 
freely, after which they are imaged by absorption jl^/ , as 
shown in Fig. 1 for k within the phonon regime. Two 
adjacent clouds are seen. The left cloud corresponds to 
the condensate. The right cloud, displaced in velocity 
space by hk/m, corresponds to the excitations. During 
the time-of-flight, the density of the atomic cloud de- 
creases dramatically, resulting in a decrease in the speed 
of sound. The phonons of wave number k are thus trans- 
formed ^ into free particles, which are subsequently im- 
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onant frequencies for the left and right directions, which 
removes the effects of the Doppler shift resulting from 
any sloshing of the condensate in the trap during the 
Bragg pulse. 




FIG. 1. The Bragg and condensate clouds, (a) Average 
of 2 absorption images after 38 msec time-of-flight, following 
a resonant Bragg pulse with fc = 2.8 fim~^ , in the phonon 
regime, (b) Cross-section of the same image. The dashed 
line is a Gaussian fit to the condensate cloud, used to find 
the zero of momentum. The radial and axial coordinates are 
indicated by r and z, respectively. 

To determine the efficiency of stimulation of excita- 
tions by the Bragg pulse, the total momentum in the ax- 
ial direction relative to the center of the condensate cloud 
is computed from the image, in the combined regions of 
the two clouds Js^ . The total momentum is divided by 
Nohk, where A^o is the average number of atoms in the 
condensate during the Bragg pulse, to obtain the effi- 
ciency. This efficiency is somewhat exaggerated though, 
because the total momentum includes momentum from 
the release process. 

The thus-measured efficiency P(fc,a;) for each direc- 
tion is shown in Fig. 2, for k — 2.8 fim^^. The curve of 
P{k,uj) is well-approximated by a Gaussian plus a con- 
stant. The constant results from the background in the 
images. The symmetric shape of P{k,uj) about the res- 
onance frequency probably reflects the spectral shape of 
the Bragg pulse, rather than the intrinsic shape which is 
expected to be asymmetric Therefore, the notation 
P{k,uj) is employed, rather than S{k,uj), whose shape is 
the intrinsic shape. 

The resonant frequency is taken as the center value of 
the Gaussian fit to P(k,uj), as shown in Fig. 2. The 
excitation energy uj{k) is taken as the average of the res- 



FIG. 2. The efiiciency P{k,Lj) for k = 2.8 /im"^ The 
open and filled circles are for left and right-travelling clouds, 
respectively. The lines are fits of a Gaussian plus a constant. 

The integral of P{k, oj) over w, equal to the integral of 
S{k,Lu), is related to S{k) by ||3| 

S{k) = 2{Trn%NjBy^ J P{k, uj)duj (6) 

where 51 — F^/ {A/S.)\/ IaIb / Isat is the two-photon Rabi 
frequency, F is the line width of the 5P3/2, F = i excited 
state, and Igat is the saturation intensity. 

Fig. 3a shows the measured excitation spectrum. The 
solid line is the Bogoliubov spectrum in the LDA (|^), 
without any fit parameters, for ^ = 1.91 kHz. The mea- 
sured spectrum is seen to agree well with the Bogoliubov 
spectrum. A linear phonon regime is seen for low fc, and 
a parabolic single-particle regime for high k. The line 
at A: = separates these two regimes. The excitations 
seen to have the smallest value of oj/k are the phonons. 
Therefore, by the Landau criterion, the superfluid veloc- 
ity Vc is bounded by w/fc for the phonons. 

The inset of Fig. 3a shows the low k region of the ex- 
citation spectrum. To extract the initial slope from the 
data, is fit to the points with k less than 3 /im^^, 
with /i taken as a fit parameter. The fit is not shown 
in the figure. The result gives the speed of sound for 
the condensate to be Ce// = 2.0 ± 0.1 mm sec~^, which 
is also the measured upper bound for v^., by the Landau 
criterion. This value is in good agreement with the theo- 
retical value of 2.01 ± 0.05 mm sec~^, which is the fc = 
slope of the Bogoliubov spectrum for /i — 1.91 ± 0.09 
kHz. The line at 2ttR~^ indicates the excitation whose 
wavelength is equal to the Thomas-Fermi radius of the 
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condensate in the axial direction. The measured exci- 
tation spectrum agrees with the LDA, even for k values 
approaching this lower limit of the region of validity. For 
the two lowest fc-values measured, the wavelength of the 
excitations are long enough that they are clearly visible, 
as density modulations in the condensate cloud. 
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FIG. 3. (a) The measured excitation spectrum tj(fc) of a 
trapped Bose-Einstein condensate. The solid line is the Bo- 
goliubov spectrum with no free parameters, in the LDA for 
/i = 1.91kHz. The dashed line is the parabolic free-particle 
spectrum. The vertical line at shows the separation be- 
tween the collective and single-particle regimes. The vertical 
Une at 2nR~^ shows the lower limit of the region of validity 
of the LDA. For most points, the error bars are not visible 
on the scale of the figure. The inset shows the linear phonon 
regime, (b) The difference between the excitation spectrum 
and the free-particle spectrum. Error bars represent Icr sta- 
tistical uncertainty. The theoretical curve is the Bogoliubov 
spectrum in the LDA for /i — 1.91 kHz, minus the free-particle 
spectrum. 

Fig. 3a also shows the parabolic spectrum for free par- 
ticles, ?ifc^/(2m). The measured excitation spectrum is 
clearly above this curve, reflecting the interaction energy 
of the condensate. To emphasize the interaction en- 
ergy, Lu{k) is shown again in Fig. 3b, after subtraction 
of the free-particle spectrum. The theoretical curve is 
the Bogoliubov spectrum in the LDA, given by (|^), mi- 



nus the parabolic spectrum for free atoms. This curve 
approaches a constant for large fc, given by the second 
term in (||). 

The closed circles in Fig. 4 are the measured static 
structure factor S{k), by (^. The values shown have 
been increased by a factor of 2.3, giving rough agree- 
ment with S{k) from Bogoliubov theory in the LDA (|^). 
Eq. (Q) is indicated by a solid line. The required factor 
of 2.3 probably reflects inaccuracies in the various values 
needed to compute fi/j. The open circles are computed 
from (H), using the measured values of uj{k) shown in 
Fig. 3a. The rough agreement between the closed and 
open circles is consistent with the relation (||), within the 
multiplicative constant applied to the closed circles. For 
the determination of S(k), it is critical that the apparent 
number of excitations is not enhanced by extra momen- 
tum obtained during the release process. Therefore, the 
number of excitations is determined by the number of 
atoms in the excitation cloud, rather than by the total 
momentum. This technique fails for the two points with 
the lowest k values, where many of the atoms do not exit 
the condensate cloud. For these points, the measured 
S{k) seen in Fig. 4 is significantly reduced. 
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FIG. 4. The filled circles are the measured static structure 
factor, multiplied by an overall constant of 2.3. Error bars 
represent la statistical uncertainty, as well as the estimated 
uncertainty in the two-photon Rabi frequency. The solid line 
is the Bogoliubov structure factor, in the LDA for fj, — 1.91 
kHz. The open circles are computed from the measured ex- 
citation spectrum of Fig. 3, and Feynman's relation (^. For 
the open circles, the error bars are not visible on the scale of 
the figure. 

For large k (short wavelength), S{k) approaches unity, 
corresponding to non-interacting, uncorrelatcd atoms. 
For long wavelengths however, S{k) counterintuitively 
approaches zero. For decreasing fc, the condensate 
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contains increasing numbers of atoms with momentum 
hk. These atoms, rather than creating additional den- 
sity fluctuations with wavelength 27r/fc, actually suppress 
such fluctuations, because the atoms are correlated in 
pairs with momenta ±hk, and opposite phase 

Since S{k) is always less than unity for the values of k 
measured here, the density fluctuations are never greater 
than in the uncorrelated case. However, S{k) is expected 
to have a peak (a roton) at a wavelength comparable to 
a, which is much shorter than the wavelengths reported 
here, but the increase of S{k) at this roton is on the order 
of a^n. For a BEC in an alkali gas, a^n ^ 10^*, so the 
roton is negligible. This is in contrast to superfluid *He, 
where S{k) has a significant roton, corresponding to a 
minimum in uj{k). 

For k > 6.8 /im^^, s-wave scattering is clearly visible. 
The measured effective scattering cross section decreases 
with decreasing fc, as predicted in |?|/. 

In conclusion, we report the first measurement of the 
excitation spectrum of a Bose-Einstein condensate, and 
the static structure factor. This spectrum displays sev- 
eral features consistent with theory. The excitation spec- 
trum consists of a linear phonon regime, as well as a 
single-particle regime. The linear regime provides an 
upper limit for the superfluid critical velocity, by the Lan- 
dau criterion. The excitation spectrum agrees quantita- 
tively with the Bogoliubov spectrum, in the local density 
approximation. This agreement continues even for long- 
wavelength excitations, close to the limit of the region of 
applicability of the approximation. The density fluctua- 
tions implied by the static structure factor agree with the 
excitation spectrum within a multiplicative constant, via 
Feynman's relation. Feynman's relation is thus verified 
within an overall constant. 
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